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Abstract 

We revisit the constraints imposed by electric dipole moments (EDMs) of nucleons and 
heavy atoms on new C P- violating sources within supersymmetric theories. We point out 
that certain two-loop renormalization group corrections induce significant mixing between 
the basis- invariant CP-odd phases. In the framework of the constrained minimal super- 
symmetric standard model (CMSSM), the CP-odd invariant related to the soft trilinear 
A-phase at the GUT scale, 9a, induces non-trivial and distinct CP-odd phases for the three 
gaugino masses at the weak scale. The latter give one-loop contributions to EDMs enhanced 
by tan/3, and can provide the dominant contribution to the electron EDM induced by 6 a- 
We perform a detailed analysis of the EDM constraints within the CMSSM, exhibiting the 
reach, in terms of sparticle spectra, which may be obtained assuming generic phases, as 
well as the limits on the CP-odd phases for some specific parameter points where detailed 
phenomenological studies are available. We also illustrate how this reach will expand with 
results from the next generation of experiments which are currently in development. 
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1 Introduction 



Electric dipole moments (EDMs) of the neutron [1] and heavy atoms and molecules [2, 3, 
4, 5, 6, 7, 8] are primary observables in probing for sources of flavor-neutral CP violation. 
The high degree of precision with which various experiments have put limits on possible 
EDMs translates into stringent constraints on a variety of extensions of the Standard Model 
at and above the electroweak scale (see, e.g. [9]). Currently, the strongest constraints on 
CP- violating parameters arise from the atomic EDMs of thallium [2] and mercury [3], and 
that of the neutron [1]: 

\d T \\ < 9 x l(T 25 ecm 

|d Hg | < 2xl0~ 28 ecm (1.1) 
\d n \ < 6xl0~ 26 ecm. 

When interpreted as a quantity induced purely by the electron EDM d e , the measurement 
of g?ti can be translated into a tight bound, \d e \ < 1.6 x 10~ 27 ecm. 

While the Standard Model CKM paradigm has received remarkable overall confirmation 
with the observation of CP- violation in the mixing and decays of neutral P-mesons, the 
motivation for anticipating new sources of CP- violation is undiminished. This is in part 
through the need to explain baryogenesis, and also through the generic appearance of new 
CP- violating phases in models of new physics introduced to stabilize the Higgs sector. In 
particular, supersymmetric models with minimal field content, i.e. the MSSM, allow for 
the presence of several CP- violating phases even in the most restrictive ansatz of flavor 
universality in the squark and slepton sectors. The null experimental EDM results pose a 
serious problem for the MSSM with superpartner masses around the (natural) electroweak 
scale. Indeed, a typical CP-violating SUSY phase of order one combined with 0(100 GeV) 
masses for the superpartners would violate experimental constraints by up to three orders 
of magnitude [10]. Moreover, ongoing and planned EDM experiments are aiming to im- 
prove the level of sensitivity to underlying CP-odd sources by several additional orders of 
magnitude. 

Over the years the SUSY CP problem has led to a number of suggestions for how this 
tuning might be alleviated: 

• small phases: The simplest solution in principle is of course just to demand that the 
diagonal phases are small (or zero). However, a framework with small phases would 
run counter to the natural Standard Model interpretation of CP- violation in the 
mixing and decays of P-mesons, which requires no analogous tuning. Of course, such 
a scenario is possible, e.g. in the framework of spontaneous CP- violation [11, 12], but 
requires a significant amount of model building, particularly in terms of constraints 
on the flavor structure, in order to suppress the flavor- diagonal phases while allowing 
the effective CKM phase to be order one. 

• flavor off-diagonal CP- violation: As a principle to motivate small diagonal phases, 
one may interpret the EDM constraints, in concert with the success of the CKM 
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paradigm, as implying that certain symmetries should ensure that all CP-violation 
arises via an interplay with flavor structures. This is an appealing viewpoint, and 
can be taken as motivation to try and divorce CP- violation from the SUSY-breaking 
sector. However, even in the case that the soft terms are real, any significant non- 
universality in these terms, i.e. misalignment between the A terms and the Yukawas, 
will generically lead to large EDMs via the large CP- violating phases in the Yukawa 
matrices (see e.g. [13]). Thus, even within this seemingly rather restrictive scenario, 
further constraints must be imposed either on the flavor sector [14] or the SUSY- 
breaking sector to lessen the impact from the EDM constraints. 

• cancellations: The MSSM contains in principle many new CP-violating phases, and 
while a lot of these are highly constrained by the flavor sector, there are still several 
new flavor-diagonal phases. Indeed, there are more phases than experimental EDM 
observables with which to bound them, and thus various partial cancellations and 
degeneracies are to be expected [15, 16, 17, 18]. However, the most recent analyses 
[19, 20, 21] using a generic set of flavor-diagonal phases and the three competitive 
EDM bounds suggests that such cancellations are difficult to achieve if at all. 

• decoupling: One can clearly reduce the constraints by lifting the scale of the super- 
partners [22]. If one only considers the leading one-loop contributions to the EDMs, 
this can be done in ways that do not significantly affect the level of tuning in the 
Higgs sector, i.e. by keeping the 3rd generation scalars light, as is also required by 
cosmological constraints [23]. However, various subleading effects then come into play 
and limit the extent to which the EDMs can be suppressed. One of our goals here 
will be to map out precisely how far EDMs can reach in terms of the superpartner 
scale. 

Taken as a whole, one may observe that these scenarios could allow for various partial 
suppressions of the EDMs but, although various models exists, within the generic supersym- 
metric framework we currently lack a compelling symmetry argument for why the EDMs 
should be suppressed by many orders of magnitude as they are in the Standard Model. 
Thus, with the current level of experimental progress, we can expect many of these scenar- 
ios to be put to the test in the coming few years, providing ample motivation for further 
theoretical consideration of the EDM sensitivities. 

To summarize the current situation, recall that on restricting attention to the framework 
of the constrained MSSM [24, 25], one reduces the number of new fundamental phases to 
two, which can be chosen to be the phases of the A and /x parameters at the GUT scale. It 
is well known that the phase of fi is particularly effective in inducing EDMs, and therefore 
for reasonably generic sparticle spectra is tightly constrained, while the phase of A can in 
principle be maximal (See e.g. [15]). For example, this may be explained by the tan/5- 
enhancement of the one-loop EDMs induced by the phase of /i, while the corresponding 
enhancement of the A-phase arises only at the two-loop level. In fact, as we shall discuss 
below, this argument is far from watertight as other contributions can also be important. 
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The present work has three main goals. The first, and more technical, goal is to update 
the calculation of the EDMs by including certain two-loop corrections in the running of the 
soft-breaking parameters from the GUT to the weak scale. At first sight, the inclusion of 
the two-loop running should not substantially change the prediction for EDMs in terms of 
the SUSY phases. However, a complex A parameter induces imaginary corrections to the 
gaugino masses, which then induce the EDMs via one-loop SUSY threshold diagrams that 
are enhanced by tan (3. This significantly modifies the EDM predictions as a function of 
the A-phase at the GUT scale and remarkably enough, in the case of even moderate tan /3, 
these new contributions can provide the dominate source for the electron EDM. 

The second goal of the paper is to present an updated and complete analysis of the 
constraints imposed by EDMs on the parameter space of the constrained MSSM (CMSSM). 
We focus on the CMSSM in order to incorporate full phenomenological analyses. However, 
for various reasons it seems that relaxing the universality assumption does not in fact 
significantly alter the conclusions [20, 21], although this issue perhaps deserves further 
attention. We will update the corresponding results of [20, 21], by incorporating the new 
Tl-EDM bound and most significantly by treating the qualitatively distinct large tan/? 
regime, utilizing the results of [26]. The third and final goal is to emphasize the significant 
improvements in sensitivity, by orders of magnitude in most cases, that are likely within a 
few years due to the current and next-generation EDM searches. 

The paper is organized as follows. In the next section we detail the new contributions 
to the EDMs induced primarily by two-loop running of the soft-breaking parameters. In 
Section 3, for completeness we give a short summary of the EDM formulae that relate 
the observables to the Wilson coefficients of the CP-odd effective Lagrangian, which are 
in turn determined by SUSY diagrams. We also outline the anticipated sensitivity of the 
next generation EDM experiments. Section 4 describes the numerical analysis within the 
CMSSM framework, and presents the limits on CP violating phases as well as the potential 
reach of EDMs in terms of the SUSY mass spectrum. Some concluding remarks appear in 
section 5. 



2 RG evolution of the CP-odd invariants at two loops 



We begin by exhibiting a simplified formula for the one-loop contributions to the electron 
EDM in terms of the CP- violating SUSY phases [19], 



d P 



16tt 2 M s 2 usy 




tan (3 sin Arg^A^m^) 



+ ||sin[Arg(M 1 M e )] 



(2.2) 



where Mi(M 2 ) are the U(1)(SU(2)) gaugino masses, all the couplings and masses are nor- 
malized at the weak scale, and we chose a SUSY parameter point with the absolute value of 
all soft breaking parameters in the selectron and gaugino sector equal to Msusy- This leads 
to a transparent overall normalization for d e . The soft breaking parameter in the Higgs 
sector, mi2, enters Eq. (2.2) via the relative phase of the two Higgs vacuum expectation 
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values. It is possible to choose m\ 2 (also denoted Bji or b in an alternative notation) and 
consequently V\ and v-i to be real. It is easy to see why the phases of fi and the Wino mass 
give larger contributions to d e than does the phase of A. For example, 



d e (0„9 M2 ) 5d ^ sm(9, + 9 M2 ) sm^ + 



d e (6 Ae ,e Ml ) ~ 2gl tan(3 X sin(^ e -Z[) " X sin(^ e ~C) ' ^ 

which implies that the sensitivity to the phase of /i is enhanced by two orders of magni- 
tude relative to the phase of A, even for moderate values of tan (3. For hadronic EDMs 
this effect is also present, although the advantage of 9 ^ is less pronounced because both 
phases can generate EDMs proportional to g\. The enhancement (2.3) provides the pri- 
mary explanation for why the phase of \i is more severely constrained than the phases of 
the A-parameters. 

Even with the restrictions of flavor universality and proportionality, the number of 
independent CP-odd phases can be rather large. Working within the CMSSM scenario, 
one assumes a common phase for the trilinear parameters at the GUT scale, and a common 
phase for the gaugino masses. This reduces the number of physical CP-odd phases to two 
and a traditional choice of basis at the GUT scale comprises 6> M and 9a- 

(9 A = Arg(A f ),9, = Aj:g(pi)y, (2.4) 
Arg(M,) = 0; 

#12 = Arg(m^ 2 ) - tuned to ensure kig(m\ 2 )Ew = 0. (2.5) 

The utility of this basis lies in the invariance of the Arg(Mj) = condition under one-loop 
renormalization group evolution, and also the invariance of 6> M . If needed, the basis (2.5) can 
be rotated to any other convenient basis in the CMSSM framework by appropriate phase 
redefinitions of the matter and gauge superfields. The choice of Arg(m^ 2 ) at the GUT scale, 
which is not invariant under RG evolution [27, 16], is such as to ensure the reality of the 
Higgs vevs. Thus, choosing different values for 9^ and 9 a requires a re- adjustment of 6> 12 at 
the GUT scale in these conventions which, although aesthetically unappealing, simplifies 
calculations by eliminating the dependence on 9\ 2 of the EDMs. 



2.1 Further contributions to d e (9 a) 

With this basis in mind, consider the relative enhancement of the (6? M + 9m 2 ) contribution 
to d e as apparent in (2.3). One observes that for moderate values of tan/5 it is sufficiently 
large to counteract a further loop suppression factor. Thus, in the absence of a /i-phase, 
one may enquire whether higher- loop corrections, induced by 9a, may become sizable in 
this way at large tan/? and compete with the standard one- loop contribution to d e (9A)- 
A number of 2-loop threshold corrections to the fermion EDMs have been considered in 
the MSSM literature, notably supersymmetric variants of the Barr-Zee diagrams [28] which 
involve the Higgs sector. However, a number of other corrections of this type exist. Indeed, 
its clear that top-stop threshold corrections on the gaugino line, as shown for d d in Fig. la, 
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Figure 1: On the left is an example of a tan (3- enhanced two- loop threshold correction to dd(0A)- If A* 
is real, the phase enters via the circled left-right stop insertion. On the right we show a more significant 
tan /3-enhanced threshold correction to d e (0A) [29]. 

may give a tan /3-enhanced contribution to fermion EDMs df{9j it ) and quark color EDMs 
(CEDMs) d q {9j it ). However, its relative size compared to the conventional one-loop term, 
ignoring hierarchies in \x and A, is of order tan f3 (m t /mi) 2 / (167r 2 ), and thus not particularly 
sizable. 

These corrections would clearly be much more important in the case of d e were it possible 
for a dependence on 9 a to manifest itself via the chargino diagram, and make full use of the 
enhancement factor of (9(10 tan/?) in (2.3). Threshold corrections of this type were noted 
by Pilaftsis [29], and the most prominent example is shown in Fig. lb. The internal stop- 
bottom loop provides a complex wavefunction correction of order £ ~ vn t A* t / (167r 2 m 2 ) 
which mixes H2 with W + . Taking the limit where the stops are heavy, the resulting 
contribution to d e scales as Ad e ~ m e tan pim(A t )/((167i 2 m{) 2 fi), which despite the loop 
supression can be comparable with, or even larger than, the 1-loop neutralino diagram. 
Moreover, such contributions are often considerably larger than the more widely known 
Barr-Zee type diagrams, and deserve to be included in more general analyses. However, in 
the restricted framework of the CMSSM, while competetive with other contributions they 
are somewhat suppressed by RG effects which generically render the ratio lm(A t )mt/m 2 
rather small at the weak scale. Indeed, we will not need to study them in detail because, 
remarkably enough, it turns out that there is an even larger class of contributions which 
are less sensitive to detailed features of the spectrum. 

In particular, one can consider a diagram where in the top-stop loop, e.g. of Fig. la, we tie 
the Higgs lines together to generate a two-loop correction to the running of the gluino mass. 
While further loop suppressed, this diagram is enhanced by a large RG log and the relative 
contribution of such a diagram compared to Fig. la is of order ln(/^GUT/Mz) (jnil m t) 2 1 (167T 2 ), 
which can easily be greater than one. In the case of quark EDMs, such corrections are still 
generally negligable. However, when such RG corrections are accounted for on the Wino 
line of the chargino diagram, the corrections to d e (6A.) as shown in Fig. 2 can be sizable. 
Despite the large loop supression, it benefits from the enhancement factor in (2.3), a large 
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Figure 2: A two- loop RG correction to ImM 2 , induced by 9a which contributes to d e (0A) via the chargino 
threshold diagram (shown in the mass-insertion approximation). Despite the loop suppression, the result 
is enhanced by tan/3, a large RG log, and a ratio of couplings, and thus can be larger than the one-loop 
neutralino contribution as discussed in the text. 



RG log, and since the insertion is divergent it is not suppressed by the generically large 
squark masses. Remarkably, within the CMSSM and for moderate to large tan (3, it seems 
that this 3-loop effect can easily be larger than the leading 1-loop threshold correction to 
d e (0A), and also the 2- loop threshold corrections in Fig. lb. Thus it will have important 
consequences for the EDMs, as we will now discuss more quantitatively below. 



2.2 Two-loop RG evolution and EDMs 

To be more precise, at two-loop order, imaginary corrections to the gaugino masses may be 
generated if ImA is nonzero. The corresponding part of the RG equation for Mi [30], 

Jt lm{Ml) = f T, de C lMA f )Tr(Y\Y f ), (2.6) 

with (in the DR scheme) 

/ 26/5 14/5 18/5 \ 
C{= 6 6 2 , (2.7) 
V 4 4 ) 

spans all the flavors, but is dominated by the top, and for large tan/3 bottom and tau, 
contributions. 

While we will not attempt a full treatment of the two-loop RG equations, we will take 
into account the relevant contributions that generate Im(Mj) at the weak scale. In order 
to estimate the order of magnitude of this effect, we can consider a simple scenario where 
we focus on the evolution of M 2 , and retain only the largest (top) Yukawa term. One then 
finds, (d / dt)lm(M2) — 31^ 2 a2/(167r 3 )Im(y4), which to leading-log order integrates to, 

lm(M 2 (M z )) ~ - 3a2 3 ( 2 ^ 3 UT) Im^o) In ~ -6 x KT 3 Im(,4 ), (2.8) 

where we have inserted the GUT-scale value for a 2 , but used the weak scale value for Y t . 
This estimate, which accounts only for the top-stop contribution, provides a reasonably 
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accurate approximation for most regimes. The two-loop suppression factor indicates why 
resummation of the large logs does not change the result significantly. Written in basis- 
independent form, the estimate (2.8) expresses corrections to three weak scale CP-odd 
invariants, Arg^Mjm 2 ^), as a function of the GUT scale invariants Arg(ilPA). 

As noted above, a particular example of the relevant contributions to the EDMs is 
shown explicitly in Fig. 2. It combines the two-loop running of the soft parameters with 
the one-loop SUSY threshold contribution to the EDM. It is clear that to a certain degree 
the effect of Im(Mj) mimics Im(/i) inside such a diagram, and thus is most prominent when 

— * 

#gut = (Ofi, 6 a) — (0, 9a)- These contributions are clearly also present for quark EDMs and 
CEDMs, but as discussed above are less competitive as the corresponding enhancement 
factor is an order of magnitude smaller than (2.3). 

The result (2.8) has a significant effect on the calculation of EDMs induced by 6a for 
several reasons: 

• The two-loop running destroys the universality of the gaugino mass phases at the 
weak scale, which therefore have to be taken into account in EDM calculations. 

• The resulting size of the imaginary parts of Mj appears to be small, O(10~ 2 — 10~ 3 ) x 
MgusY (2-8). However, in combination with the enhancement factor ~ 10 tan (3 (2.3), 
the resulting correction to d e can be 0(1 — 10), for moderate-to-large tan/3, as com- 
pared to the "standard" result for d e (6 A ). 

• The two-loop RG induced correction to d e (6A) has the opposite sign relative to the 
standard result, providing a regime where the two contributions can cancel each other. 

Putting together the various numerical factors, from RG-running and the SUSY threshold- 
induced EDM, we see that the new contributions induced by Im(Mj(6U)) correspond to a 
three-loop correction which is enhanced by 0(1 — 10) tan f3 x ln(//Q UT /M§), and as em- 
phasized above is often larger than the leading 1-loop threshold correction and various 
additional 2-loop threshold corrections. The question then arises as to whether there are 
other three-loop effects with a comparable enhancement, i.e. oc tan 2 /3. It is well known, 
however, that such contributions do not exist. Enhancement by a single power of tan/5 
results from d e oc y e (H 2 ) oc m e tan (5. Any additional power of tan (5 can arise only from ad- 
ditional insertions of the Yukawa couplings of charged leptons or down-type quarks, which 
cannot lead to a numerical enhancement. 



3 An EDM summary 

In this section, we briefly recall the dependence of the observable EDMs on the most 
significant flavor-diagonal CP-odd operators at 1 GeV (see [31] for a recent review). Up to 
dimension six, the corresponding effective Lagrangian takes the form, 

q 2 - 

£ eff = -M*- 6 G a llv G^ a 
n 32tt 2 ^ 
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i=e,u,d,s i=u,d,s 

+\w f ahc G%G^ b G p ^ + £ C^m^l^) + ■■■ (3-9) 

For the four-fermion operators, the sum runs over all light fermions, and one could in prin- 
ciple include additional tensor structures [32, 33, 34, 35]. In practice, since these operators 
require a double helicity flip, implying additional {m q /vEw) 2 suppression factors, they are 
generally negligible. The exception to this rule is at large tan/5, where grow as (tan/?) 3 , 
and they can compete with the other sources in certain parameter regimes [36, 26]. 

In supersymmetric models with CP-violating phases larger than ~ 10~ 8 in the flavor- 
diagonal sector the only reasonable strategy to avoid the strong CP problem is to postulate 
the existence of the Peccei-Quinn relaxation mechanism [37] , which removes the theta term 
from the list of contributing operators. We will adopt this strategy here. 

We will now turn to a quick synopsis of the dependence of the observable EDMs on the 
remaining dimension five and six operators. The physical observables can be conveniently 
separated into three main categories, depending on the physical mechanisms via which an 
EDM can be generated: EDMs of paramagnetic atoms and molecules, EDMs of diamagnetic 
atoms, and the neutron EDM. 



3.1 EDMs of paramagnetic atoms — thallium EDM 

Among various paramagnetic systems, the EDM of the thallium atom currently provides 
the best constraints. Atomic calculations [38, 39, 40] (see also Ref. [9] for a more complete 
list) link the atomic EDM with d e and various CP-odd electron-nucleon interactions, of 
which we shall only consider the most relevant, namely Csei^eNN, 

d Tl = -5854 - e 43 GeVCp 9 , (3.10) 

where we have furthermore retained only the leading isosinglet component of Cs- The 
relevant atomic matrix elements are calculated to within 10 — 20% accuracy [41]. 1 

The dependence of Cs on the four fermion sources C e , for i — d,s, b, can be expressed 
in the form [36, 26], 

m9 _ 29 MeV k x 220 MeV 66 MeV(l - 0.25k) 

— ^de r Ls se hC/fe , l-J- iJ J 

m d m s m b 

where k = {N\m s ss\N) /220 MeV, and we will adopt the value obtained using NLO x?T, 
k ~ 0.50±0.25 [42]. Larger values for the pion-nucleon sigma term a n N, advocated recently 
(see e.g. [43]), would result in a somewhat larger value for k ~ 1.55. 

1 It is worth emphasizing that, due to the contributions from Cs, a bound on more than one paramagnetic 
species is required to infer a direct bound on d e . Thus, a conservative and model-independent bound on d e 
cannot be derived from the bound on cZti alone, and would be ~ two orders of magnitude weaker than the 
commonly quoted value obtained by dropping the dependence on Cs in (3.10). 
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3.2 Neutron EDM 



The calculation of the neutron EDM in terms of the Wilson coefficients of Eq. (3.9) rep- 
resents a difficult non-perturbative problem, and thus the precision currently obtainable 
is not competitive with the results outlined above for d^\. Nonetheless, d n plays a crucial 
role in constraining CP-odd sources in the quark sector, and we will make use of the re- 
sults obtained using QCD sum rule techniques [44, 45] (see [46, 47, 19, 48] for alternative 
approaches using chiral techniques), that we briefly recall below. 

In the presence of the PQ mechanism, it is natural to expect that the dominant con- 
tribution to the neutron EDM comes from the EDMs and CEDMs of the light quarks. 
Within the sum-rules framework, PQ relaxation also suppresses the sea-quark contribution 
at leading order, and leads to the following result [44] : 

d n (d q , d q ) = (1.4 ± 0.6)(d d - 0.25<4) + (1.1 ± 0.5)e(d d + 0.5d u ). (3.12) 

The quark vacuum condensate, (qq) = (225 MeV) 3 , has been used in this relation - the 
proportionality to d q (qq) ~ iTi q (qq) ~ fw 771 ^ removes any sensitivity to the poorly known 
absolute value of the light quark masses. Here d q and d q are to be normalized at the hadronic 
scale which we assume^ to be 1 GeV. Note also that the quark masses used for the SUSY 
calculations of d q and d q should be taken at the weak scale, where their numerical values are 
smaller than the low energy values by a factor of ~ 0.35, e.g. m d (M z ) ^ 9.5 MeV x 0.35. 

The dimension six sources, w and C qiq2 , also contribute to the neutron EDM but are 
more problematic to handle within the sum-rules framework. One can nonetheless obtain 
estimates that are useful for assessing the regimes in which the contributions in (3.12) are 
dominant. One finds [45], 

d n {w) ~ 20 MeV x e w. (3.13) 

This estimate is assessed to be valid within a factor of 2-3 [45]. When tan (3 is large, 
the dominant contributions to w arise via threshold corrections from the CEDMs of heavy 
quarks [49, 50], e.g. 5w(d b ) ~ g^d b (mb)/(327r 2 mb). 

In general there are also many possible contributions from four-fermion sources. How- 
ever, since these contribute only at large tan/?, the most relevant are Cbd — —Cdb and 
C s b — —Cbs- In this case, one can obtain an estimate of the form [26], 

dJCbd) ~ e 0.7 x 10~ 3 GeV 2 — , (3.14) 

which again is primarily useful as a means to estimate regimes where large corrections to 
(3.12) are possible. 

3.3 EDMs of diamagnetic atoms — mercury EDM 

EDMs of diamagnetic atoms, i.e. atoms with total electron angular momentum equal to 
zero, also provide an important test of CP violation [9]. The current limit on the EDM of 
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mercury [3] furnishes one of the most sensitive constraints on SUSY CP- violating phases 
[19]. However, the calculation of dn g is undoubtedly the most difficult as it requires QCD, 
nuclear, and also atomic input. 

The atomic EDM of mercury arises from several important sources, namely, the Schiff 
moment S [51], the electron EDM d e , and also electron-nucleon interactions (see, e.g. Ref. 
[9] for a comprehensive review). Schematically, the mercury EDM can be represented as 

(d h C qiq2 )\, C s [C qe ], C P [C eq ], d e ) , (3.15) 

where g n NN collectively denotes the CP-odd pion-nucleon couplings, and Cs and Cp denote 
respectively the couplings eij^eNN and eeNi^^N . The Weinberg operator does not provide 
any appreciable contribution to du g because its contribution to g n NN is suppressed by an 
additional chiral factor of m q /l GeV ~ 10~ 2 . 

A number of the atomic and nuclear components of the calculation have been updated 
recently, specifically du g (S) [52] and S(g w NN) [53]. We will also make use of an updated 
result for g n NN(di) obtained by combining chiral and sum-rules techniques [54], 

£L(d q ) ~ 4l«(4 - ^[GeV- 1 ], (3.16) 

where the relatively large overall uncertainty arises from significant cancellations between 
the leading contributions. Combining this result with the contribution of four-fermion 
operators [26], and the atomic and nuclear parts of the calculation, we obtain 

d Hg = 7 x 10~ 3 e(d u - d d ) + 1(T 2 4 

-1.4 x KT 5 eGeV 2 (®Ji£*L + 3 .3/^ + ^(1 - 0.25k)) (3.17) 

V 777.(2 m s m b / 

+3.5 x l(T 3 GeVeCs + 4 x 10~ 4 GeVeC P . 

The contributions in the second and third lines are significant in the MSSM only for large 
tan p. 

The most valuable feature of d Hg is its sensitivity to the triplet combination of CEDM 
operators di, which surpasses the sensitivity of the neutron EDM to this combination by a 
factor of a few. Moreover, despite the fact that the overall coefficient is only known with 
relatively poor precision, the dominant dependence on the (d u — dd) combination ensures 
that these uncertainties enter as an overall factor and therefore do not significantly alter 
the shape of the unconstrained band of the parameter space in the 6* M — 9a plane. 

3.4 Future experimental sensitivity 

The experimental situation is currently very active, and a number of new EDM experiments 
are already in development which promise to significantly improve the level of sensitivity 
by several orders of magnitude in the next few years. As part of our analysis in the next 
section, we will present the reach of these next-generation experiments within the CMSSM, 
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given the following assumed sensitivity in the three classes of EDMs on which we have been 
focussing: 



\Ad e \ < 3xl(r 29 ecm 
Ad D \ < 2xl(T 27 ecm 
|AdJ < 1 x 10~ 27 ecm. 



(3.18) 



These sensitivities were chosen conservatively according to expectations for the first physics 
runs of the relevant experiments discussed below. It is possible that the ultimate sensitivity 
will be be at least an order of magnitude better in most cases. 

The first constraint shown in (3.18) corresponds to a relatively conservative estimate of 
the sensitivity achievable within three active projects utilizing respectively two polarizable 
paramagnetic molecules, PbO [8] and YbF [7], and also a solid state system [55]. The 
first two of these primarily obtain their enhanced sensitivity from the fact that the internal 
electric field, to which the electron EDM is sensitive, is further enhanced by the polarization 
of the molecule. One should bear in mind that, as for Tl, these experiments do not directly 
bound d e and electron- nucleon operators may also be relevant. PbO and YbF are polarized 
by the external field in a nonlinear way and thus do not have EDMs as such. However, the 
induced shift in the Larmor precession frequency can still be written in the form, 



We are not aware of any calculations for the C^-dependence in this case, while the quoted 
sensitivity in (3.18) comes from ignoring this term and using estimates for the effective 
field S e s [56] and the sensitivity to uj^ [8, 7]. As alluded to above, if no further systematic 
issues arise, the final sensitivity could in fact be several orders of magnitude better than 
this quoted result [8, 55]. Given that the relevant atomic calculations are only partially 
completed, we will simply scale up the result for Tl and restrict ourselves to a relatively 
low value of tan/5, in order to ameliorate the problem of unknown corrections from Cs- 
Nonetheless, one should be aware that the molecular matrix element relating Cs and u^, 
although enhanced by the same factor Z 3 will not be a simple rescaling of the Tl result, and 
thus the final dependency uj^(d e , Cs) may be somewhat different from d^i(d e , Cs)- Clearly 
further progress on atomic calculations for this system would be welcome. 

The second constraint in (3.18), refers to the proposed search for the deuteron EDM 
using the BNL storage ring [57]. This proposal seems our best chance for significant im- 
provement in the diamagnetic sector. The deuteron is not a diamagnetic system of course, 
but the EDM turns out to be primarily sensitive to g^NN (assuming 6 is removed by PQ 
rotation) [58, 59, 60], 



and thus it plays a similar role in terms of its sensitivity to underlying CP-odd sources. 
Indeed, it has the significant advantage in this regard of not requiring complicated many- 
body nuclear calculations. 



hAu L ~£ eS d e + 0(C s ). 



(3.19) 



d D ~d n + d p - (1.3 ± 0.3)e^ 7V ((i g )[GeV- 1 ] 



(3.20) 
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The final constraint in (3.18) refers to the expected (initial) sensitivity of several exper- 
iments currently in development, at LANSCE [61], at ILL [62], and at PSI [63], searching 
for an EDM of the neutron. The ultimate sensitivity could again be an order of magnitude 
better. 

4 Numerical Results and EDM constraints 

Before we turn to the numerical analysis of EDMs within the CMSSM framework, for 
completeness we now list the relevant SUSY contributions to the operators in (3.9) that 
must be included: 

4.1 SUSY contributions to CP-odd higher-dimensional operators 

Within the effective theory formulated at the weak scale, or rather the mass scale of the 
superpartners, the CP-odd operators in (3.9) can be induced by the following threshold 
corrections, all of which we will include in our analysis. 

Fermion EDMs and CEDMs 

• 1-loop: At one-loop order the electron EDM receives contributions from x~ — v e and 
X° — e~ loops (see e.g. [17]). The EDMs and CEDMs of quarks receive analogous 
contributions, and in addition have generically dominant contributions from gluino- 
squark loops (see e.g. [17]). As we argued in Section 2, these one loop results must 
be formulated in the basis where ImM, at the weak scale are taken into account. 

• 2-loop: At 2-loop order, non-negligible Barr-Zee-type Higgs-mediated graphs con- 
tribute to EDMs and CEDMs [28]. The induced photon couplings, e.g. AFF and 
HFF, arise through sfermion and chargino loops, with dominant contributions from 
stop (and, for large tan/3, sbottom and stau [26]) loops. These diagrams can be 
significant when the first generation sfermions are heavy. 2 

Weinberg operator 

• 2-loop: The Weinberg operator [64] also contributes to d n via two-loop stop-top- 
gluino contributions (see e.g. [17]), particularly in the regime where the first gener- 
ation sfermions are heavy. As noted earlier, the Weinberg operator may also receive 
non-negligible corrections from heavy quark CEDMs under RG evolution to low scales 
[49, 50, 26]. 

2 Note that additional 2-loop threshold corrections, such as Fig. lb [29], can in principle be larger than 
the Barr-Zee type contributions, and thus a full analysis of 2-loop threshold corrections appears warranted. 
However, within the CMSSM, we have verified that the complex wavefunction-induced correction £ [29] 
to the chargino mass matrix is, even when non- negligible compared to other conventional sources, still 
generally subleading relative to complex 2-loop RG corrections to the 1-loop diagrams. 
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Four-fermion operators 

• 1-loop: We account for the leading tan /^-enhanced contributions [36], as computed 
in [26], induced by scalar or pseudoscalar Higgs exchange, with a tan /5-enhanced 
coupling of down-type fermions to H 2 induced by various 1-loop threshold corrections. 
There are also similar Higgs exchange contributions with loop-induced H — A mixing 
that are generically sub-dominant. 

These contributions are quite generic in any SUSY model of new physics at the TeV 
scale. We turn now to a numerical analysis of the EDM sensitivity to phases and SUSY 
masses within the framework of the CMSSM. 

4.2 EDM reach and constraints on CP-odd phases in CMSSM 

In this section we analyze all three observables, dxi, d n and du g , within the CMSSM, defined 
by the following set of universal SUSY parameters at the GUT scale: 



where mi/ 2 and m are the GUT scale universal gaugino and scalar masses respectively. The 
magnitude of the /^-parameter and the pseudoscalar mass are determined by the radiative 
electroweak symmetry breaking conditions. We will present the analysis separately for each 
of the two physical phases, 9 ^ and 9 a, which are present under these assumptions. 

Let us recall that RG running from the GUT scale in this scenario introduces consid- 
erable mass splittings in the spectrum of superpartners; for example, the gluino becomes 
much heavier than (roughly triple) the rest of the gauginos, while the squarks also become 
quite heavy, m^ q (Mz) — itlq + 6m^ 2 + 0(M§). Of particular importance here, as empha- 
sized above, is that the physical phases also run, and mix, under the RG. The two-loop RGE 
for the \x parameter is real, and therefore the phase of ji does not run; its low energy value 
is equal to its input value at the GUT scale. On the other hand, the phases of all of the 
A parameters do run and are primarily affected by the 3rd generation Yukawa couplings. 
As discussed earlier, although we are free to set the phases of the gaugino masses to zero 
at the GUT scale, the 2-loop RGEs induce phases for M\^^ when the phases of A are not 
zero. The mixing of physical phases that this induces will be illustrated below. 

In what follows, we present our numerical results for the EDMs in the context of the 
CMSSM. 

4.2.1 Phase constraints in benchmark scenarios 

Detailed phenomenological studies have found it useful to concentrate on very specific 
CMSSM parameter choices. To this end, several sets of benchmark points have been estab- 
lished [65, 66, 67]. A set of proposed post-LEP benchmark scenarios [65] were chosen to 
span the CMSSM parameter space in regions where all phenomenological constraints are 




(4.21) 
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satisfied. These include obtaining a cosmological relic density with the limits established 
by WMAP [68]. Several points were chosen to lie in the 'bulk' region at small mi/ 2 and 
mo, while others are spread along either the \ ~ t coannihilation 'tail' at larger mi/ 2 or in 
rapid-annihilation 'funnels' where the relic density is controlled by s-channel annihilation 
of the LSP through the heavy Higgs scalar and pseudo-scalar at large m 1( / 2 and m at large 
tan/5. 

We begin our analysis, therefore, with a discussion of EDMs in the context of four 
selected benchmark scenarios. Point B is a bulk point with (mi/ 2 , m , tan /3) = (250, 75, 10). 
The value of m is shifted slightly (+15 GeV) from [67] due the larger value of the top 
quark mass (we adopt m t = 178 GeV as opposed to 175 GeV used previously) and to 
several improvements in the spectrum code, most notably the inclusion of the full set of 
2-loop RGEs. Note that here we have also chosen A = 300 GeV as opposed to 0, to allow 
for effects with non-zero phases. This also adds to the shift in m . In Fig. 3a, we show 
the magnitude of the EDMs for benchmark point B in a plane defined by the phases of 
A and /i. Within the blue shaded region the Tl EDM is less than or equal to its current 
experimental limit. The blue (dashed) contours show the values of phases when the Tl EDM 
exceeds its experimental bound by a factor of 10 and 50. The black (dashed) curve within 
the shaded region has g?ti — 0. Similarly, within the red shaded region the neutron EDM is 
within its experimental bound. The red dotted curve corresponds to a neutron EDM which 
is 10 times the experimental limit. Finally, the green shaded region displays the Hg EDM 
bound and the green (solid) curves shows where the Hg EDM is 5 times its experimental 
limit. As expected, the bounds on 6^ are far stronger than those on 6 a- However, because 
the EDMs cut the phase planes differently, when limits to all three EDMs are applied, a 
relatively strong limit {\9a/k\ < 0.08) can be obtained in this case. The limit on \9^\ is 
approximately 0.0027T. 

In Fig 3b, we show the result for benchmark point D, which is a point along the x ~ t 
co-annihilation tail and (in the absence of phases) has ji < with (m^, m , tan/3) = 
(525, 130, 10), applying a 20 GeV shift in m . In this case, the sparticle masses are generally 
larger than found for point B, primarily due to the increase in mi/ 2 . Here, we do not show 
the shaded region obtained from the Hg EDM as it would fill the displayed plane. The Hg 
contours in the upper right and lower left of the panel show where the EDM is equal to its 
experimental limit. For point D, we find no limit to 6a, while — ir\ < 0.0657T. 

Point L (shown in Fig. 3c) is also a co-annihilation point, at large tan/9. Point L is 
defined by (mi/ 2 , m , tan (5) = (450, 355, 50). In this case, because tan f3 is large, we applied 
a larger shift (55 GeV) in m to obtain the correct relic density when phases are absent. 
Once again, we recover a bound on 6a, albeit a weak one of about 0.2-7T. In this case, we have 
a very strong constraint on 6^, \6^\ < 0.00057T. Finally, point M, shown in Fig. 3d, is also a 
large tan/3 point found high in the funnel region with (mi/ 2 , mo, tan/3) = (1500, 1100,57). 
Because the funnel region is very sensitive to tan/3, mi/ 2 and m t , we increased tan/3 from 
50 to 57, and lowered mi/ 2 from 1840 GeV to 1500 GeV. As for point D, we do not show 
the shaded region for the Hg EDM (as it would fill the plane as displayed). As one can 
see, there is no bound on 6a, and despite the large values of m 1 / 2 and m , there is still a 
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Figure 3: The Tl (blue dashed), neutron (red dotted), and Hg (green solid) EDMs relative to then- 
respective experimental limits in the 9^, 9 a plane for benchmark points B, D, L and M. Inside the shaded 
regions, the EDMs are less than or equal to their experimental bounds. Each of the EDMs vanish along 
the black contour within the shaded region. 
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benchmark point (m^, mo, tan/3) 


\0aM 


\0JA 


B (250,75,10) 


< 0.08 


< 0.002 


D (525,130,10) 




< 0.07 


L (450,355,50) 


< 0.2 


< 0.0005 


M (1500,1100,57) 




< 0.004 



Table 1: A summary of the constraints on the phases within the four benchmark scenarios shown in 
Fig. 3. All points use \A \ = 300 GeV, and the quoted bound on the phase of /j, for point D is modulo tt. 

significant bound on |^| of 0.0047T. 

The constraints for all four benchmark scenarios are summarized in Table 1. 

4.2.2 General constraints on the CMSSM parameter space 

We now turn to a more general analysis of the CMSSM (m^, mo) plane for fixed values of 
tan/3. We will present the results in several different ways. One motivation for restricting 
our attention to the CMSSM is that we can subject the theory to a full suite of different 
phenomenological constraints. On the following parameter space plots we will indicate the 
excluded and preferred regions due to the following: 

• LSP dark matter: Regions in parameter space where the LSP is charged, generically 
a stau, are not viable and can be excluded. Moreover, if one demands that the relic 
density is consistent with the range Vlh 2 = 0.0945 — 0.1287 determined by WMAP 
[68], this limits the viable region to thin strips in the (mi/2,wio) plane. Furthermore, 
the relic density is generically too large on one side of the WMAP strip. We will 
calculate this relic density in the absence of phases, but will restrict the magnitude 
of the phases so that their impact should be negligible [15, 16]. 

• Contributions to b — > S7: Within the flavor sector, constraints on the supersym- 
metric contributions to b — > 57 are often important for low values of my 2, particularly 
at large tan f3 [69] . 

• LEP chargino and Higgs searches: The constraints resulting from the null LEP 
results for Higgs [70] and chargino [71] searches also impinge on the CMSSM param- 
eter space for small mi/ 2 . 

Our motivation is not to search specifically for regions allowed by all of these constraints, 
but to emphasize the complementarity of the constraints, and the relative importance of 
those imposed by EDMs. The resulting conclusions should then be relevant, qualitatively 
at least, in more general SUSY-breaking scenarios. 

In presenting our results, we will consider the contribution of the two phases separately, 
and we will also use two phase values of 7r/20 and it/ 6, and show contour plots of d/d exp , 
the EDM relative to its current experimental bound (or future level of sensitivity). These 
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phases are both within the linear regime and thus the contours can easily be converted to 
bounds on the phases, 9 < 9 / (d/d exp ), where 9 denotes the nonzero phase chosen for the 
plot. This allows for a simple rescaling of our results given updated future limits. 

The smaller phase value, 7r/20, was chosen to restrict attention the phenomenologically 
most interesting region, where the spectrum is not highly tuned. Note that the tuning of 
this phase is comparable to the minimal tuning that must be imposed on the spectrum in 
the CMSSM, as implied by the LEP bounds and the electroweak precision tests (see e.g. 
[72] for a recent analysis within the MSSM). We have also included results with a "near 
maximal" phase, tt/6, to illustrate the full reach of the EDM constraints with regard to 
the parameter space of the CMSSM and the SUSY spectrum. The value n/6 was adopted 
as the boundary of the linear regime, to allow for a simple interpretation of the contours, 
and, if viewed as a "maximal" phase, also takes into account a conservative (factor of two) 
uncertainty in some of the hadronic EDMs. 

We will now present the results in more detail. 
Low tan (3 and future sensitivity 

In Fig. 4a,b, we show the current constraints in the CMSSM plane for tan f3 — 10 due 
to the Tl, neutron, and Hg EDMs for 9 A = 0.05tt (a) and 9^ = 0.05tt (b). As in Fig. 3, 
the blue (dashed), red (dotted), and green (solid) curves correspond to the Tl, neutron, 
and Hg EDMs respectively. As described above, the contour labels show the magnitude of 
the EDM relative to the current experimental bound. Throughout the (m^jMo) plane, 
we have taken m t = 178 GeV, m b (m b ) = 4.25 GeV, and \A \ = 300 GeV, with 9 A = 0.05vr, 
9 p, = in a) , and 9 ^ = 0.057T, 6 A — in b). In each of the panels, the brown shaded region 
in the lower right corresponds to that portion of the CMSSM plane where the right-handed 
stau is the LSP, and is therefore excluded. The green shaded region in the left of each 
panel is also excluded as there the supersymmetric contribution to b — > 57 is too large. 
Also shown are the current constraints from LEP based on the masses of the chargino 
(nearly vertical black dashed line) and Higgs boson (red dot-dashed) . Regions to the left of 
these curves are excluded. The grey shaded wisp-like region, corresponds to that portion 
of the plane where the cosmological relic density lies in the range Qh 2 = 0.0945 — 0.1287 as 
determined by WMAP. Any area above this shaded region is excluded as the relic density 
is too large. Note that the very thin area below the WMAP strip and above the stau LSP 
region is allowed so long as there is another source for dark matter. 

In Fig. 4a, we show the contours for which each of the three EDMs considered, saturates 
its experimental bound (curves labeled 1). Regions to the left of (or below) these curves 
are excluded. The constraints from d n and c?H g are relatively uniform. However, we see that 
for the Tl EDM, there is a line through the plane where a cancellation occurs, and d^i = 0. 
Interestingly enough, despite the relatively low value for tan f3, this cancellation line can be 
understood through the impact of 2-loop RG corrections to ImM,. 

To see this, recall that for low tan/3 the Tl EDM is dominated by 1-loop threshold 
corrections to d e , which to a reasonable approximation take the form shown in Eq. (2.2). 
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a) tan B = 10 , 9 A = 0.05tz, A n = 300 GeV W tan B = 10 , 0„ = 0.05ti, A = 300 GeV 
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c) tan B = 10 , A = 0.05ti, A n = 300 GeV d) tan B = 10 , 6„ = 0.05ji, A = 300 GeV 
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Figure 4: The Tl (blue dashed), neutron (red dotted), and Hg (green solid) EDMs relative to their current 
experimental limits (a,b) and the electron (blue dashed) , the neutron (red dotted) , and the deuteron (black 
solid) EDMs relative to possible future experimental limits (c,d) in the toi/2,«io plane of the CMSSM. In 
each panel, |Ao| = 300 GeV with tan/3 = 10. The shaded regions are described in the text. In a) and c), 
6 A = 0.05tt, and in b) and d), 8 M = 0.05tt. 
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Since /i is real, the negative RG corrections to ImM 2 can lead to a cancellation between the 
two leading contributions. From (2.8), 



d e oc 



5^tan/3//Im(M 2 ) + 2gflm(M*A e ) , (4.22) 



using the leading-log estimate for ImM 2 in (2.8), and ignoring the subleading imaginary 
correction to Mi, we can estimate the condition for dominance of the first term in (4.22) as 

where in the second relation we have used ~ {~ m2 H 2 ) 1 ^ ~ l-6^i/2, M ~ Oim^, 
and lm(A e ) ~ 1.6Im(A ) to provide an indication of the size of tan/5 needed to achieve a 
cancellation in the relevant part of the (mi/ 2 ,m ) plane. We see that moderate values of 
tan/3 are already sufficient, and indeed this cancellation is apparent in Fig.4a. Its precise 
shape depends on a number of other contributions that we have neglected, which become 
more important at large tan/? as we will discuss below. For completeness, the plot also 
indicates contours for the EDMs which are 5 times their experimental bound, and one 
contour for Tl where the EDM is 10 times it experimental value (the innermost of the blue 
dashed curves). 

In contrast, when 9 ^ = 0.057T, as in Fig. 4b, the constraints become much stronger. All 
the curves are now labelled by the value of d/d exp , i.e. the ratio by which the predicted 
EDM exceeds the experimental bound. For this choice of 9^, we see that the entire region of 
the plane shown in the Figure is excluded. The Tl EDM, which now exhibits no cancellation 
line as it is dominated by the first term in (2.2), supplies by far the strongest constraint. 
Nonetheless, d n and dn g also provide strong constraints across the plane. As described 
above, since the phase is small and the contribution to the EDM is thus linear, we can also 
use the curves to place a limit on 9^. For example, the region to the right of the Tl curve 
labeled 50 is allowed so long as 9^ < 0.057r/50 = 0.001-7T. Similarly for the other curves in 
this panel. 

In the lower half of Fig. 4, we consider how the present situation will be modified 
with the next generation of EDM experiments, and exhibit the anticipated sensitivities. 
In Fig. 4c and d, we show d e (blue dashed), the neutron (red dotted), and the deuteron 
(black solid) relative to the anticipated future level of sensitivity as described in section 3.4. 
We emphasize that there are several experiments in this list sensitive to d e , and we have 
simply deduced the sensitivity to d e from these, ignoring corrections from electron-nucleon 
interactions. For small values of tan j3 as we have here, this should not be a significant 
problem but one should bear in mind that these conclusions cannot easily be extrapolated 
to other regimes, e.g. to large tan (3. 

In panel c), the contours shown are 5, 10, and 50, where the latter are unlabeled. 
For d e , contours of and 1 are also shown. As one can see we expect the deuteron and 
neutron to give comparable limits which are quite complementary to that obtained from 
the future bound on the electron EDM from PbO, YbF, etc. As is apparent from the 
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comparison of Fig. 4a and Fig. 4c, the limits on the SUSY parameter space are expected 
to be dramatically improved by future experiments, particularly so for the dependence on 
9 a- For 9a = 0.057T, the entire (displayed) plane would be excluded by both the n and D 
EDMs, provided that future experiments will see zero EDMs, thus forcing a smaller value 
for the phase of A . In panel d), we see that the d e limits become so strong that even at 
the endpoint of the co-annihilation region at mi/ 2 ^ 900 GeV, we are forced to very small 
phases, 6> M < 5 x 10~ 5 7r. 

Constraints for large tan (3 

The large tan f3 regime exhibits a number of new qualitative features as concerns the de- 
pendence on the underlying sources. The 1-loop ^-dependence of the down-type fermions 
is linearly enhanced, while as we have discussed RG corrections to ImMj can similarly 
provide significant corrections for the dependence on ^. A number of 2- loop Barr-Zee 
type contributions [28, 73] to the EDMs and CEDMs are also tan /3-enhanced and become 
significant in this regime, while the largest qualitative change arises through the fact that 
various four-fermion operators provide sizable contributions [74, 36, 26], and compete with 
the constituent EDMs. The most prominent example of this is the contribution of Cs to 
g?ti as described in Section 3. A comprehensive study of the relative importance of these 
contributions was carried out in [26], and here we will present combined results for the 
EDMs within the framework of the CMSSM. It is apparent from the results of [26] that the 
restrictive ansatz of the CMSSM does not exhibit all the possible regimes that may arise 
in a more general MSSM scenario, and new qualitative features are restricted primarily to 
the behaviour of d^\. 

Our results are presented in the upper two rows of Fig. 5, for tan (3 = 30 and 57 respec- 
tively, and again for cases with either nonzero 9a = 0.057T or 9 ^ = 0.057T. As mentioned, the 
hadronic EDMs, d n and dn g , are still primarily dominated by the one-loop contributions 
and thus are not altered significantly for nonzero 9a (cf. Fig. 5a,c with Fig. 4a), while 
they scale up linearly with tan/5 for nonzero 9^ (cf. Fig. 5b, d with Fig. 4b). We recall that 
although the 2-loop RG-induced corrections for nonzero 9a are present and tan /3-enhanced, 
they are not particularly large for the quark EDMs and CEDMs as all contributions arise 
predominantly via squark-gluino loops. Indeed, it is apparent from Fig. 5a and c that, for 
our rather small value of 9a = 0.05-7T, the EDMs are less constraining than the limits from 
b — > S7. On the plots for tan f3 = 57, one may note the appearance of the funnel-like region 
where dark matter annihilations are mediated by s-channel heavy Higgs exchange. 

Let us now focus attention on the contours for d^. At first sight, the results for tan (3 = 
30 may not look dramatically different from those for tan/3 = 10. We see once again a 
cancellation line extending out now to somewhat larger values of m^. However, its precise 
form is now determined by a more complex interplay between many different effects. Indeed, 
from (4.23) we might expect that the ImM2-induced corrections would now dominate the 1- 
loop contribution to d e . This is indeed the case, but for tan (3 = 30 the 2-loop corrections and 
the contributions from Cs are also competitive, and the observed cancellation line results 
from the interplay of all of these sources. It is only for sufficiently large tan (3, as observed 
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a) tan P = 30, 9 A = 0.05n, A = 300 GeV b) tan P = 30, 6 = O.OStc, A = 300 GeV 
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Figure 5: As in Fig. 4, the Tl (blue dashed), neutron (red dotted), and Hg (green solid) EDMs relative 
to their current respective experimental limits in the toi/2,too plane of the CMSSM for tan/3 = 30 (a, b, 
and e), and tan/? = 57 (c, d, and f). In each panel, |^4o| = 300 GeV. In a, c, e, and f), 6a = 0.057T, and 
in b) and d), 0^ = 0.057T. In panels e and f) only 1-loop RGEs were used to compute the supersymmctric 
particle spectrum. 
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for tan (3 = 57 in Fig. 5c, that the cancellation line disappears, as the various contributions 
fall out of balance, although there are still some partial intermediate cancellations [26] . 

We have not included plots outlining the future sensitivity in this case, due to the lack of 
atomic calculations detailing the contributions from four-fermion operators for the relevant 
paramagnetic sources. Nonetheless, by comparison with Fig. 4c it should be clear that 
much of the displayed plane will again be covered by the EDM reach, even for 9a = 0.057T, 
providing a significant improvement over the current results in Fig. 5a and c. 

Impact of two-loop RG-evolution of the phases 

To gauge the effect of the inclusion of the 2-loop RGEs, we have also shown in Fig. 5e 
and f the constraints on the CMSSM plane when only the 1-loop RGEs are run for non- 
zero 6a- For clarity, we have only displayed contours with d = d exp . As anticipated in our 
discussion in Section 2, the Tl EDM is most affected by the inclusion of the 2-loop RGEs 
and as one sees, the cancellation line at tan/? = 30, which using the 1-loop RGEs barely 
appears on the top left corner of the plot, is rotated clockwise and shifted to close to the 
middle of the plane. Similarly, for tan (5 = 57, the cancellation that is apparent when using 
the 1-loop RGEs is shifted into the stau LSP region in Fig. 5c. 

One should bear in mind in viewing these plots that the spectrum is also slightly per- 
turbed in shifting from 1-loop to 2-loop RGEs, but this is a minor effect and the dominant 
change in the EDMs is due to the induced corrections to ImMj. The significant shift in the 
WMAP band is due to the fact that the funnel region is highly sensitive to the pseudoscalar 
Higgs mass. 

Full EDM reach in the CMSSM parameter space 

In the analysis above, we have chosen rather small phases, consistent with the level 
of implied tuning in the CMSSM, to focus in on the phenomenologically most interesting 
region. However, it is also clearly of interest to determine precisely how far the EDMs can 
reach in terms of the parameter space of the CMSSM, and more generally in terms of the 
sparticle spectra. Thus, we will now boost the phases to n/6 ~ 0.177T, which we will take as 
the boundary of the linear domain, so that we can still use the figures to approximate the 
limit on the given phase at any specific point in the CMSSM plane. As discussed earlier, 
if 7r/6 is taken as a maximal phase value, it also provides a conservative treatment of some 
of the calculational uncertainties. 

Our results are shown in Fig. 6, where the current (a and b) and future (c and d) 
contours are shown for both 6a (a and c) and 9^ (b and d). Here we show only the 
value of the EDM relative to its experimental limit and suppress the other cosmological 
and phenomenological constraints. Where possible, contours of 0, 1, 5, 10, 50, and 100 
are shown. Note that the shaded region in the upper left corner of each panel above the 
solid black curve is not viable in the CMSSM as there are no solutions which allow for 
radiative electroweak symmetry breaking. (The boundary zone near this line is often called 
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Figure 6: As in Fig. 4, the Tl (blue dashed), neutron (red dotted), and Hg (green solid) EDMs relative 
to their current experimental limits (a,b) and the electron (blue dashed), the neutron (red dotted), and the 
deuteron (black solid) EDMs relative to possible future experimental limits (c,d) in the m 1 / 2 ,mo plane of 
the CMSSM for tan/3 = 10. In each panel, \A \ = 300 GeV. In a and c), 9 A = 0.17tt, and in b) and d), 

Of, = 0.177T. 
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the focus point region, or the hyperbolic branch.) When 9a = 0.17n and 9^ = 0, we 
see from panel a), that the reach is not terribly strong. For comparison, we show the 
expected reach of the LHC based on its ability to detect colored sparticles with masses 
up to 2 TeV [75, 76]. However, we observe that the future sensitivity of the neutron and 
deuteron EDM experiments will extend beyond the reach of the LHC for m < 4.5 TeV. 
While it is interesting that the cancellation line apparently weakens the constraints from 
paramagnetic sources, we emphasize that corrections from four-fermion sources may change 
this significantly in the case of PbO and YbF, and that this cancellation may not survive 
in more general variants of the MSSM. 

In contrast, when 9 ^ = 0.17n, the current reach of all three EDM experiments already 
exceeds the projected reach of the LHC as seen in panel b) of Fig. 6. For completeness, we 
show in panel d) the expected reach from future EDM experiments. This figure shows that 
for large 9^ the EDMs can be sensitive to superpartner mass scales as high as a few tens of 
TeV. While such scenarios are apparently highly tuned, one may try to rescue naturalness 
via SUSY breaking mechanisms, outside the CMSSM framework, where only the first two 
generation scalars are pushed above 10 TeV in order to cure the SUSY flavor problem by 
'brute force', while keeping the 3rd generation scalars light in order not to exacerbate the 
fine-tuning problem in the Higgs sector [22]. It is clear that the next generation EDM 
searches will directly probe such scenarios. 

The shape of these contours is easily understood given the large RG-induced admixture 
of mi/2 to the squark masses. We see that the contours for the hadronic EDMs are squashed 
towards the m axis as a consequence of the rapid suppression of the EDMs as my 2 becomes 
large. This effect is absent for the paramagnetic EDM bounds, and we observe a more 
uniform shape for these contours. 

5 Concluding Remarks 

We have performed a general analysis of the EDM constraints on the CMSSM parameter 
space, combining a number of weak-scale threshold contributions to the fermion EDMs and 
CEDMs, and also four-fermion operators. The latter contributions are generally subleading, 
but become important at large tan j3 and must be taken into account for a complete analysis 
of this regime. We have emphasized two main points as part of this analysis. Firstly, as the 
experimental situation is rapidly evolving, we have presented estimates for the sensitivity 
of the next generation of EDM experiments in each of the three sectors which currently 
provide the strongest constraints; namely paramagnetic sources sensitive to the electron 
EDM, and the neutron and deuteron which are sensitive to (different) combinations of 
quark/gluon operators. The most striking impact of the expected levels of improvement in 
sensitivity is that the reach of EDMs induced by 9 A increases significantly, to the extent 
that it exceeds that of the LHC for large phases. The reach in terms of 9^ already probes 
superpartner mass scales that cannot be directly accessed at the LHC. 

The second point, of technical importance, that we have emphasized is that 2-loop RG- 



24 



induced corrections to ImMj can significantly enhance the dependence of the electron EDM 
to 9a, particularly so at large tan/3. While on first inspection it would appear that these 
RG-induced contributions generically dominate the Tl EDM by a factor of a few, it turns 
out that the relative signs are such that significant cancellations occur, and the bounds 
are often weakened, at least for moderate tan/5. For large tan/3, a number of different 
contributions become competitive, but on closer inspection one may observe by comparing 
Fig. 5c with f that the constraints from Tl do in fact become stronger with the inclusion of 
the RG-induced corrections, although overall they remain relatively weak in this regime. A 
summary of some of the constraints on the phases in specific benchmark scenarios is shown 
in Table 1. 

A natural question is whether these conclusions will hold in more general examples of 
the MSSM, where a number of the universality assumptions are lifted, and the number of 
independent phases is correspondingly increased. While the most recent explorations of 
this kind [20, 21, 26] suggest that the picture is not dramatically modified, this is clearly a 
question worthy of further study. 
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